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Abstract: The reconstruction of extensive air showers (arrival direction, core position and 
energy estimation) by the surface detector of the Pierre Auger Observatory is discussed to- 
gether with the corresponding accuracy. We determine the angular reconstruction accuracy 
as a function of the station multiplicity by using two different aproaches. We discuss statisti- 
cal and systematic uncertainties in the determination of the signal at 1000 m from the core, 
S(IOOO), which is used to estimate the primary energy. 
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Introduction 

The Pierre Auger Observatory consists of two 
independent components: the fluorescence de- 
tector (FD) and the surface detector (SD) 
[1]. We have determined the angular reso- 
lution of events recorded by the surface de- 
tector alone, on an event by event basis, 
from the zenith (9) and azimuth (0) uncer- 
tainties obtained from the geometrical recon- 
struction, using the relation described in [2]: 
F{r]) = 1/2 {V[9] + sin^ie) V[(t>]) ,wherer] 
is the space- angle, and V[9] and V[(j)] are the 
variance of 9 and </> respectively. We define 
the angular resolution (AR) as the angular ra- 
dius that would contain 68% of showers com- 
ing from a point source, AR = 1.5 yjF{r]). 
We checked the angular resolution using the re- 
dundant information given by a sub-array com- 
posed by adjacent detectors. 
The parameter used to infer the energy of the 
surface detector events (S'(IOOO)) is studied 
and its systematic and statistical errors are de- 
termined. The event-by-event error estimation 
is checked with full Monte Carlo simulations . 
The unavoidable fluctuations in this parameter 
caused by fluctuations in the shower develop- 
ment is evaluated with simulations for different 
primary assumptions. 



Angular Resolution 

The arrival direction of a SD event is deter- 
mined by fitting the arrival time of the first 
particle in each station to a shower front model. 
The precision achieved in the arrival direction 
depends, on the clock precision of the detector 
and on the fluctuations in the flrst particle ar- 
rival time. In [3] an empirical model has been 
developed to determine the uncertainty in the 
time measurement of each individual detector 
participating in the event. 

The model of the shower front used in the min- 
imization procedure, be it spherical, parabolic, 
or even planar also influences the uncertainty 
in the arrival direction determination, but not 
as much as the time measurement precision. It 
has been shown in [3] that a parabolic model 
for the shower front adequately describes the 
data. 

On a event by event basis 

Given the two inputs: a model for the time 
variance and a model for the shower front, 
the angular resolution can be calculated on 
an event by event basis out of a minimiza- 
tion procedure. In Fig. 1, we show our an- 
gular resolution as a function of the zenith an- 
gle for various station multiplicities (circles: 3 
stations, squares: 4 stations, up triangles: 5 
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Figure 1: Angular resolution {AR) for the SD 
as a function of the zenith angle (0). The AR 
is plotted for various station multiplicities. 

stations, down triangles: 6 stations or more). 
The data used to build this plot spans from 
January/2004 to March/2007. 

As it can be seen, the angular resolution is bet- 
ter than 2° in the worst case of vertical showers 
with only 3 stations hit. This value improves 
significantly for 4 or 5 stations^. For 6 or more 
stations, which corresponds to events with en- 
ergies above 10 EeV, the angular resolution is 
in all cases better than about 1°. 

Using station pairs 

A new sub-array of pairs has been recently de- 
ployed as a part of the Surface Detector array. 
These are adjacent detectors located 11 m 
apart, and therefore are sampling the same re- 
gion of the shower front. To do this analysis, 
events with at least three pairs are selected. 
The reconstruction is then performed twice, 
each time using the time information of one 
of the tanks in each pair. This provide two 
quasi-independent estimates of the geometry. 
In Fig. 2 we show the space-angle difference 
between these two estimates for showers with 
3, 4, and 5 or more stations. 

The distributions are then fitted to the ad- 
justed Gaussian resolution function {dp oc 
g-i? /2<T c?(cos(77)) dr], where 77 is the angle 
between the two reconstructions of the same 
shower) to obtain a. The angular resolution 
(68% contour), which is given by 1.5 times cr, is 
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Figure 2: Space-angle difference between two 
SD estimates of the event geometry for differ- 
ent multiplicities (see text for more details). 

in agreement with the one obtained on a event 
by event basis. 



Energy Estimator 

The surface detector only samples the proper- 
ties of an air shower at a limited number of 
points at different distances from the shower 
axis (r). An observable has to be then defined 
to estimate the shower size. To avoid the large 
fluctuations in the signal integrated over all 
distances caused by fluctuations in the shower 
development, Hillas [4] proposed to use the sig- 
nal at a given distance {S{r)) to classify the 
size of the shower. In Fig. 3 we show the pre- 
dictions from Monte Carlo simulations of the 
magnitude of the fluctuations in S{r = 1000) 
as a function of zenith angle. The relative fluc- 
tuations are found to be independent of energy 
and its magnitude is ~ 10% for most of the 
cases studied. 

The experimental error in the estimation of the 
signal size at a given core distance depends on 
the spacing of the array. In [5] it has been 
shown that for the Auger array spacing the 



1. For 4 and 5 stations the AR is very similar 
because in the fitting procedure they have the same 
number of degrees of freedom. 
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Figure 3: Relative spread due to shower to 
shower fluctuations for different compositions 
(blue-iron, red-proton, black-mixed composi- 
tion). 

optimum distance {ropt) to minimize this ex- 
perimental error is ~1000 m. Therefore, the 
observable that we use to relate to the pri- 
mary energy will be the signal size at 1000 m 
(5'(1000))2. However, it should be noted that 
Topt fluctuates from event to event and in- 
creases to larger core distances (~ 1500 m) 
when there are saturated stations [5]. 

To estimate S'(IOOO) it is necessary to adopt 
a lateral distribution function (LDF) that de- 
scribes the fall-off of the signal size with the 
distance to the shower axis. The function used 
here is a modified NKG function given by: 
S{r) = 5(1000) (j^)-^ (T^r^ where r is 
the distance to the shower axis in meters, S{r) 
is the signal size at a core distance r, 5(1000) 
is the size parameter of the shower, and /3 is 
called the slope of the LDF. 

S(IOOO) uncertainties 

The signal sizes in each station are then used 
to estimate the core location and 5(1000), with 
(3 being a fixed parameter. The fitting error 
in 5(1000) is a consequence of the uncertainty 
of the observed signal size largely due to the 
finite dimension of the detectors. This will be 
the statistical error in 5(1000) (o-s(fooo))- The 
uncertainty in the signal sizes has been mea- 
sured directly using pairs of stations located 
close to each other in the field [6] . 

The second source of error in 5(1000) is a 
systematic (ct^^^qqq^ arising from the lack of 
knowledge of the true LDF shape for a par- 
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Figure 4: The average systematic and sta- 
tistical error in S(IOOO) as a function of 
log 5(1000). The data has been divided in two 
sets (events with- without stations saturated). 

ticular event. If the rgpt of a given event is 
close to 1000 m, the fitted 5(1000) is indepen- 
dent of the value of f3 assumed [5] . When it is 
not, fluctuations in the event by event /3 give 
rise to a systematic error. The value of (3 to 
be used in the reconstruction has been esti- 
mated empirically: in a small subset of events 
(5(1000) > 20 VEM and having more than 5 
stations) the (3 is left as a free parameter as 
well. We then parameterize the fitted values of 
/3 as a function of zenith angle and 5(1000). 
The deviation from this parameterization is 
calculated for each event and the RMS (ap) 
parameterized as a function of 5(1000) (no 
dependence on zenith angle has been found). 
The resuh is the following: (7^3(5(1000)) = 
0.71 X exp(-0.976 log(5(1000))). We then re- 
peat N times the reconstruction of each event, 
fixing (3 to values sampled from a Gaussian dis- 
tribution centered around the prediction with 
the sigma given above. The RMS of the fit- 
ted 5(1000) from these set of fits is then the 
systematic error of 5(1000) (f 5^1000) ) ■ 
In Fig. 4 we show the average systematic and 
statistical error of 5(1000) as a function of 
log(5(1000)). The data has been divided in 
two sets: events with (without) saturated sta- 



2. S'(IOOO) is measured in units of VEM, i.e. the 
signal produced by a vertical centered muon. 
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Figure 5: Total error in 5(1000) calculated on 
an event by event basis from the data. The 
data is separated in two sets: events with 
(without) saturated stations. The lines cor- 
respond to the predictions from full MC calcu- 
lations (see text for details) 



tions. Two features are clearly seen: a) the 
error in events with saturated stations is sys- 
tematically 5% larger, b) crs{iooo) dominates 
the error budget for 5(1000) <40 VEM. No de- 



pendence of Cg^^pQp-j 
has been found. 
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Using Full Monte Carlo Simulations 

To benchmark our error estimation wc have 
created a library of Corsika showers for pro- 
ton primaries, zenith angles 9— 0-12-25-36- 
45-60 degrees and energies logj^Q £'(ey)=17.8- 
18.0-18.2-18.4-18.6-19.0-19.5-20.0. For each 
Corsika shower, we calculate the true S'(IOOO) 
and it is then used to generate 10 (25) events 
(depending on the energy) with random core 
positions. 

The reconstruction procedure used for the data 
is then applied to the simulations. For each 
zenith angle and energy we fit the distribu- 
tion of log ^ 5'(io°o°o°lriie ^ to a Gaussian func- 
tion. The mean value and sigma are then pa- 
rameterized as a function of 5(1000)*''"'^. No 
zenith angle dependence has been found. A 
bias in the reconstructed 5(1000) is only found 
for 5(1000) < 10 VEM. The sigma of this dis- 
tribution is the quadrature combination of the 



statistical and systematic error in 5(1000). In 
Fig. 5 we show the comparison of the sigma of 
these distributions with the average total error 
obtained on an event by event basis. The data 
is separated in two sets: events with (without) 
saturated stations. The circles correspond to 
the total error obtained on a event by event 
basis, the lines are the predictions from full 
Monte Carlo simulations. The agreement is 
excellent except for a slight overestimation of 
the error (~4%) for saturated events at large 
energies. 

Conclusions 

The angular resolution of the surface detec- 
tor was determined experimentally, checked us- 
ing the pairs data set and found to be bet- 
ter than 2° for 3-fold events {E < 4 EeV), 
better than 1.2° for 4- folds and 5-folds events 
(3 < < 10 EeV) and better than 0.9° for 
higher multiplicity events {E > 10 EeV). 

The error of the parameter used to infer the en- 
ergy of the surface detector events (5(1000)) 
has been determined experimentally, checked 
using full Monte Carlo simulations and found 
to be better than 8% (12%) at the highest en- 
ergies for events with (without) saturated sta- 
tions. At high energies, the fluctuations in 
5(1000) are dominated by fluctuations in the 
shower development. 
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